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STABILITY  CONSIDERATIONS  FOR  LIGHT-ION  BEAN 
TRANSPORT  IN  Z-DISCHARGE  CHANNELS 


1.  INTRODUCTION 

The  development  of  terawatt-level  ion  beams  has  generated  a  great  deal 
of  interest  in  using  light-ion  beams  to  drive  thermonuclear  pellets. 1_^ 

Target  design  studies  for  light-ion  beams  indicate  that  ^  2  MJ  of  ions  must 
be  delivered  in  t  10  nsec  to  an  ^  1  cm  diameter  pellet  in  order  to  achieve 
high-gain  thermonuclear  ignition.^"0  Since  present  technology  can  provide  up 
to  5  TW  single-generator  modules  from  which  up  to  200  kJ  of  ions  can  be 
extracted  in  'v  50  nsec,  a  multimodule  system  is  required.  In  addition,  a 
transport  scheme  and  a  method  for  beam  pulse  compression  are  needed. 

One  possible  transport  scheme  involves  the  use  of  a  z-discharge  channel 
for  transporting  a  prefocused  ion  beam  (Fig.  I).11  Focusing  is  achieved  by 
a  combination  of  geometric  and  magnetic-field  focusing  prior  to  injection  into 
the  channel.12'13  Beam  pulse  compression  results  from  ramping  the  diode 
voltage  such  that  the  tail  of  the  beam  catches  up  to  the  front  of  the 
beam.13'14  The  ideal  diode  voltage  waveform  is  <J>(t)  =  4>0 (l-t/ta) "2  for 

0  £  t  <  tb  <  ta  where  tb  is  the  beam  pulse  duration  and  ta  is  the  beam  arrival 
time  at  the  target. 

Assuming  that  the  z-discharge  channel  is  produced  in  a  MHD  stable  con¬ 
figuration,  the  question  of  the  effects  that  the  passage  of  the  beam  will  have 
on  the  equilibrium  and  stability  of  the  beam-plasma  system  is  an  important 
one.  The  MHD  response  of  the  plasma  has  been  treated  elsewhere13  and  will 
only  be  briefly  reviewed  here.  Analysis  of  stable  beam  propagation  in 
straight  and  tapered  channels,  as  well  as  in  bumpy  channels  (subject  to 
sausage  instability)  has  also  been  done  previously.18  This  work  shows  that, 
in  the  absence  of  microinstabilities  driven  by  the  beam,  good  beam  transport 
and  bunching  is  possible  under  the  conditions  set  by  MHD  considerations 
(which  will  be  outlined  in  Sec.  2).  However,  in  the  presence  of  microinsta¬ 
bilities,  beam  transport  and  bunching  can  be  seriously  affected. 

Analysis  of  electrostatic  (ES)17  and  electromagnetic  (EM)18  velocity- 
space  instabilities,  which  can  grow  on  a  time  scale  much  faster  than  the  beam 
pulse  duration,  will  be  reviewed  in  this  report.  The  problem  will  be  con¬ 
sidered  with  the  goal  of  identifying  the  conditions  for  good  transport  and 
bunching.  The  results  will  be  presented  in  a  general  form  so  that  they  may 
be  applied  to  beams  of  various  low  atomic  number  species  propagating  in 
channel  plasmas  of  different  compositions. 


Manuscript  submitted  on  May  13,  1981. 
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2.  BEAM-PLASMA  SYSTEM  IN  THE  TRANSPORT  CHANNEL 

The  beam-plasma  system  consists  of  a  focused  ion  beam  propagating  down 
the  axis  of  an  externally-driven  z-discharge  plasma  channel . , 19  The  ion 
beam  is  focused  at  the  entrance  to  the  plasma  channel  (see  Figs.  1  and  2)  with 
velocity  components  transverse  to  z  given  by  V^/Vz  =  tan  a  <<  1.  The  current 
being  driven  in  the  preformed  z-discharge  channel  provides  the  radial  confine¬ 
ment  of  the  beam.  Because  of  the  small  perpendicular  beam  energy,  the  channel 
current  can  be  much  less  than  the  beam  current.  A  high  plasma  density  in  the 
channel  insures  good  beam  charge  neutralization . 20  Good  beam  current 
neutralization  in  the  interior  of  the  beam  also  occurs,  so  that  the  total 
magnetic  field  is  comparable  to  that  associated  with  the  preformed  channel 
established  before  beam  injection.  Because  Jp  =  Jch  -  Jb  *  -Jb»  the  electron 
drift  velocity  is  approximated  by  Ve  -  nbzbvb/^pnp<  where  n^  and  np  are  the 
beam  and  plasma  ion  densities,  and  Zp  are  the  beam  and  plasma  ion  charge 
states,  and  Jp,  J^  and  Jch  are  the  plasma,  beam  and  preformed  channel  current 
densities,  respectively.  Note  that  the  electron  density  is  Zpnp. 

Hydrodynamic  modeling  of  the  background  plasma^  shows  that  a  uniform 
net-current  model  is  appropriate  for  the  early  times  associated  with  passage 
of  the  beam  front.  This  is  because  the  low-temperature  channel  is  established 
microseconds  before  beam  injection  so  that  complete  magnetic  diffusion  occurs. 
Later  in  the  ion  pulse,  Jp  x  B  expansion  of  the  beam-heated  high-temperature 
plasma  reduces  the  magnetic  field  strength  in  the  interior  of  the  channel. 

The  built-up  field  in  the  expanding  cylindrical  shock  wave  is  also  enhanced 
by  significant  current  non-neutralization  in  the  cool  plasma  surrounding  the 
beam-heated  channel.  The  maximum  field  strength  just  outside  the  ion-beam 
radius  can  exceed  that  established  by  the  preformed  z-discharge  current  by  a 
large  factor.  Thus,  at  late  times  during  beam  passage,  the  magnetic  field 
distribution  can  be  approximated  by  a  surface-current  model.  Although  the 
induced  Vr  x  b@  electric  field  is  important  when  considering  beam  energy 
losses  during  transport,  at  no  time  does  the  electric  field  become  large 
enough  to  significantly  affect  the  stability  analysis. 

The  linearized  stability  analysis  presented  here  strictly  applies  only 
when  the  mode  under  consideration  grows  on  a  time  scile  faster  than  any 
changes  in  the  beam-plasma  system  (i.e.,  y  >  t^1)-  For  9rowth  on  a  time 
scale  slower  than  the  beam  pulse  duration,  the  perturbation  analysis  breaks 
down  since  small  perturbations  will  be  washed  out  by  the  zero-order  changes 
in  the  system  before  growth  can  occur.  Thus  for  those  modes  which  are  shown 
to  e-fold  less  than  once  during  the  passage  of  the  beam,  the  results  should 
be  interpreted  as  showing  that  significant  growth  does  not  occur. 

The  distribution  of  particles  in  axial  velocity  is  illustrated  in  Fig.  3. 
The  plasma  ions  form  a  stationary  background  while  the  drift  of  the  plasma 
electrons  provides  for  beam  current  neutralization.  Because  of  the  high 
plasma  density,  this  drift  velocity,  Ve,  does  not  exceed  the  electron  thermal 
velocity  even  before  the  beam  heats  the  plasma.  As  the  beam  passes  through  a 
given  point,  the  plasma  is  heated  and  the  conductivity  increases.  The 
electron-ion  collision  frequency,  ve^,  decreases  but  the  electrons  generally 
remain  collisional  (ve^  >  uuce)  at  all  times  during  the  pulse  for  purposes  of 
the  stability  analysis.  Since  the  electron-ion  equilibration  time  is  on  the 
order  of  the  beam  pulse  duration  the  electron  temperature  will  not  exceed 
the  ion  temperature  by  more  than  a  factor  of  ten. 
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The  ion  beam  density  is  typically  much  less  than  the  plasma  density  and, 
beams  of  interest  are  nonrelativistic .  The  spread  in  axial  velocities, 
designated  by  Vt  in  Fig.  3,  is  on  the  order  of  Vfc,aj/2  before  beam  bunching 
occurs.  Here  is  the  beam  velocity  and  c^,  is  the  maximum  ion  injection 
angle  into  the  channel.  For  am  =  0.1-0. 2  rad,  V^/V^  is  relatively  small.  As 
the  beam  bunches,  V^/Vb  car  increase  by  an  order  of  magnitude  as  the  faster 
ions  generated  later  in  the  pulse  catch  up  to  slower  ions  at  the  front  of  the 
beam.  The  beam  pulse  duration,  Tb,  decreases  and  the  beam  density,  nb, 
increases  as  the  beam  bunches.  The  dependence  of  the  beam  parameters  on  axial 
position  due  to  bunching  will  be  expressed  explicitly  as  Vt(z),  x^Cz),  etc. 

Since  the  ES  stability  analysis  depends  strongly  on  the  shape  of  the 
beam  distribution  function,  it  is  important  to  use  a  theoretical  mouel  which 
contains  the  appropriate  physics.  For  a  distribution  in  vz  which  has  a 
large  slope  on  the  low  velocity  side  such  as  illustrated  in  Fig.  4a,  a 
resistive  instability  persists  in  spite  of  the  thermal  spread  in  vz  and  the 
high  frequency  of  collisions  between  the  plasma  electrons  and  ions.2^»22 
However,  if  the  slope  of  the  distribution  function  is  not  as  sharp,  such  as 
for  the  Gaussian  distribution  shown  in  Fig.  4b,  electron-ion  collision  are 
damping  if  the  thermal  spread  in  vz  is  sufficiently  large.  The  slope  of  the 
distribution  function  on  the  high  velocity  side  does  not  affect  the  stability 
analysis. 


Before  the  beam  bunches,  the  distribution  in  vz  is  determined  by  the 
injection  condition.  Since  the  ions  are  injected  into  the  channel  nearly 
uniformly  over  a  range  of  angles  predominately  in  the  r-z  plane  and  since 
the  time-averaged  axial  velocity  for  an  ion  injected  at  a  given  a  varies  like 
Vb(l-a^/2),  the  distribution  in  vz  rises  slowly  as  indicated  in  Fig.  4c. 
Additional  smoothing  out  of  the  distribution  in  vz  will  result  from  beam 
energy  spreading  due  to  radial  variations  in  the  diode  voltage.  The  stability 
properties  of  this  distribution  closely  resemble  those  of  the  Gaussian  distri¬ 
bution.  The  beam  ion  distribution  after  bunching  is  also  similar  in  shape 
to  Fig.  4c.  Since  the  ion  diodes  which  are  used  in  the  experiments  have  a 
constant  impedance  behavior  during  the  duration  of  the  ion  pulse,  more  ions 
(higher  currents)  are  generated  at  higher  energies  as  the  diode  voltage  ramps 
upward.  Thus  in  the  bunched  state  there  are  more  ions  at  higher  velocities 
than  at  lower  velocities  and  the  spread  in  axial  velocities  is  on  the  order 
of  (2Zde /m^)  .  Here  Z<j  is  the  charge  state  of  the  beam  ions  in  the 

diode  and  4>0  and  4>f  are  the  initial  and  final  diode  voltages.  This  spread 
can  be  considerably  larger  than  the  initial  spread  before  bunching  (V^VyX^ /2) . 


Thus  a.  Gaussian  distribution  can  be  used  to  properly  model  the  ion  beam 
distribution  function  both  before  and  after  bunching  when  considering  ES 
modes.  When  the  analysis  is  not  sensitive  to  the  detailed  shape  of  the  dis¬ 
tribution  function  (e.g.  for  analyzing  the  stability  of  EM  modes),  simpler 
models  may  be  used. 


3.  LINEAR  STABILITY  ANALYSIS 

17  18  22 

The  results  of  linear  stability  analyses  reported  previously  ’  * 

are  presented  here  in  a  more  general  form  so  that  they  may  be  applied  to 
beams  of  various  species.  The  composition  of  the  channel  plasma  is  also 
unspecified.  Growth  rates  will  vary  as  the  beam  heats  the  channel  plasma 
and  as  the  channel  expands.  The  changes  in  the  hydrodynamic  structures  of 
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the  channel  occur  gradually  over  the  pulse  duration  of  the  beam.  Growth  rate 
expressions  appropriate  near  both  the  front  and  the  tail  of  the  beam  are 
presented  in  order  to  determine  how  hydrodynamic  changes  alter  the  results. 
Most  notably  these  changes  include  (1)  a  reduction  in  the  e-i  collision 
frequency  at  the  tail  of  the  beam  where  the  plasma  is  heated,  and  (2)  a  reduc¬ 
tion  in  the  magnetic  field  in  the  interior  of  the  channel  at  late  times  in  the 
pulse  as  the  now  highly  conducting  heated  plasma  expands  due  to  J  x  B  forces. 
Near  the  front  of  the  beam,  the  beam  ions  follow  betatron-like  orbits  while 
at  the  tail  of  the  beam,  the  beam  ions  move  in  more  straight  line-like  orbits 
with  reflections  off  of  the  magnetic  field  piled  up  at  the  edge  of  the 
channel. 


3.1.  ES  Modes 


The  ES  modes  are  the  fastest  growing  modes  when  driven  unstable. 

Unstable  growth  will  generally  reach  nonlinear  saturation  on  a  time  scale  much 
faster  than  the  beam  pulse  duration.  These  modes  involve  either  charge 
bunching  (e-b  and  e-i  modes)  or  density  bunching  (ion-acoustic  mode)  and 
generally  exhibit  their  largest  growth  for  k  =  kez,  which  is  a  result  of  the 
small  thermal  spread  in  the  axial  direction  before  axial  beam  compression 
occurs. 


3.1a.  e-b  Two  Stream,  Mode 

The  e-b  two  stream  instability  is  driven  by  the  relative  streaming 
between  the  beam  ions  and  plasma  electrons.  In  general  the  beam  can  be 
considered  warm  because  Vt(z)/Vb  >  2[mezbnb (2) /mbzpnp]1//3  even  before  beam 
bunching  occurs.  Here  mb  is  the  beam  ion  mass.  Near  the  front  of  the  beam 
the  betatron  motion  of  the  ions  can  reduce  the  growth  rate  of  the  mode  which 
is  given  by 


where  w2b  =  4Tre2Zj?)nb  (z) /mb  and  w|,e  =  4tre2Zpnp/me  are  the  beam  and  electron 
plasma  frequencies  respectively.  The  electron-ion  collision  frequency  is 
given  by  vei  =  1. 45*10-6ZpnpAei/Tj/2 ,  where  Te  is  in  eV  and  is  the 
Coulomb  logarithm.  The  reduction  factor,  R,  equals  the  fraction  of  beam  ions 
which  can  effectively  participate  in  the  wave  growth.  If  Yj  >  2ui^,  which  is 
the  usual  case,  all  ions  can  participate  and  R  =  1.  Here  wg  =  (vbwcb/rc) 15  is 
the  ion  betatron  frequency  and  wcb  is  the  beam  cyclotron  frequency.  If 
Yj  5  2wg,  only  those  ions  with  Az  <  1/k  can  participate,  where  Az  is  the 
amplitude  of  the  betatron  oscillations  about  z  =  vzt  (see  Fig.  5) .  The  ^ 
amplitude,  Az,  depends  on  the  ion  injection  conditions,  r(z=0)  and  a(z=0). 

If  this  amplitude  is  large,  the  ion  moves  across  many  wave  fronts  before  the 
instability  e-folds  even  once,  and  thus  cannot  effectively  participate  in 
wave  growth.  The  value  of  R  must  be  calculated  from  the  actual  distribution 
function.  For  a  Gaussian  distribution  R  ^  8wgVb/kVb  for  k  >  SupV^/Vt  and 
R  =  1  for  k  <  8wgVb/Vb.  At  the  tail  of  the  beam,  betatron  effects  are  less 
important,  so  that  R  =  1  in  Eq.  (1). 
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Beam  heating  of  the  plasma  can  considerably  reduce  the  damping  term 
in  Eq.  (1)  as  e-i  collisions  become  less  frequent.  Beam  bunching,  on  the 
other  hand,  reduces  the  driving  term  in  Eq.  (1)  as  vt  increases  dramatically. 
Thus  the  potentially  most  dangerous  position  for  wave  growth  lies  at  the  tail 
of  the  beam  (R=l)  at  the  beginning  of  the  transport  channel  before  signifi¬ 
cant  beam  bunching  occurs. 


3.1b.  e-i  Two  Stream  Mode 


The  e-i  two  stream  mode  is  driven  by  the  relative  streaming  between 
the  background  electrons  and  ions.  The  electrons  drift  with  an  average 
velocity  Ve~nbZbVb/Zpnp  relative  to  the  stationary  ions  and  in  general  Ve  <  ue. 
Here  ue  is  the  thermal  velocity  of  the  electrons.  The  expression  for  the 
growth  (damping)  rate  is 


where  cs  =  (T^ir^)  is  the  ion  sound  speed,  XDe  =  (Te/4Tre2Zpnri) ’  is  the 
electron  Debeye  length  and  where  a  simple  Krook  model  was  usecl  for  the 
collision  term.  In  deriving  Eq.  (2),  it  was  assumed  that  ve^  <  kcs/(l+k2x£e) **. 
When  Te  'v-  T^  this  reduces  to  ” 
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Since  Ve  is  usually  less  than  Uq,  the  mode  is  typically  stable. 

3.1c.  Ion-Acoustic  Mode 

If  Te  >>  T^,  then  the  drifting  electrons  can  drive  an  ion-acoustic 
instability.  In  this  case  Eq.  (2)  reduces  to 
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This  predicts  instability  (y  >  0)  ,  if  V  /cs  >  (l+k2x£e)-^  and  if  vei  is 
sufficiently  small.  However,  the  severity  of  the  condition  on  Te/Ti  is  often 
overlooked.  For  Ve/cs  as  large  as  5,  Te/Ti  must  be  greater  than  12  for 
instability  even  for  kXDe  <<  1.  Figure  6  shows  the  critical  values  of  Te/T^ 
for  instability  versus  k  for  various  values  of  Ve/cs  (assuming  ve^  =  0) . 

Since  Tg/T^  is  not  expected  to  reach  such  high  levels,  in  general  ion- 
acoustic  turbulence  is  not  expected. 
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3. Id.  ES  Stability  Conditions 

Since  the  e-b  mode  is  the  only  ES  instability  which  could  be  gen¬ 
erated,  it  is  important  to  state  under  which  conditions  it  may  be  avoided. 
For  stability  (Yj  <  0 )  Eq.  (1)  states  that 


z5/2n3/2  x  , 

3.2  X  107 

t3/2 
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zbVz)v 

U  ' V2(z) 
b  t 


2 

b 


(5) 


is  required.  Here  y^,  =  mjj/mjj,  Te  is  measured  in  eV  and  all  other  variables 
are  in  Gaussian  units.  This  condition  used  in  conjunction  with  the  condition 


np  *  1.4  x  1012  ^-2  lIb(z)Tb(z)/s))3  ,  (6) 

derived  from  MHD  considerations15,  set  constraints  on  n^  and  Te  for  good 
transport.  Here  Icyj  and  I)-,  are  measured  in  amps  and  S  is  the  beam  stopping 
power  of  the  plasma  measured  in  erg-cm2/g.  Eq.  (6)  gives  the  ion 
density  in  the  channel  required  for  minimum  beam  energy  loss  in  the 
channel  during  transport.  Given  this  required  ion  density,  Eq.  (5)  sets  an 
upper  limit  on  Te,  beyond  which  electrostatic  turbulence  will  set  in.  For 
good  transport  both  of  these  conditions  must  be  satisfied.  Generally  the 
stability  criterion  in  Eq.  (5) ,  is  most  severely  tested  at  the  tail  of  the 
beam  where  Te  is  largest  and  at  the  beginning  of  the  transport  channel  before 
significant  bunching  occurs  where  Vt(z)  is  smallest  (Vt(z)/v^  'v  a^/2) . 

3.2  EM  Modes 

The  EM  modes  are  slower  growing  modes  than  the  ES  modes,  so  that  even 
if  the  mode  is  unstable,  growth  may  not  have  sufficient  time  to  reach  non¬ 
linear  saturation  during  the  pulse  duration  of  the  beam.  The  modes  are  nearly 
purely  growing  (i.e.  do  not  convect  with  the  beam)  and  involve  current 
bunching  or  f ilamentation  of  the  beam  and/or  plasma  channel.  The  wave  vector, 
k,  is  oriented  perpendicular  to  the  direction  of  beam  propagation  with 
k  =  ker  for  the  Weibel  mode  (radial  current  bunching)  and  k  =  keg  for  the 
whistler  mode  (azimuthal  current  bunching).  In  addition,  the  cylindrical 
geometry  dictates  that  kr  >  27r/rj3  and  kg  i  Jl/rj-,  where  l  is  an  integer. 

3.2a.  Beam-Whistler  Mode 


The  beam-whistler  instability  is  driven  unstable  by  the  relative 
streaming  motion  between  the  beam  ions  and  the  channel  plasma.  Since  the 
wave  vector  is  in  the  azimuthal  direction  and  the  beam  ions  execute  their 
betatron  orbits  in  the  r-z  plane,  the  betatron  motion  of  the  beam  has  little 
effect  on  the  mode.  However,  a  small  spread  in  angular  momentum,  which  is 
observed  experimentally,  can  reduce  the  growth  rate  significantly.  With  no 
spread  in  angular  momentum  (V@/c2  <  2vei“pbvb/ccope^  the  9rowth  rate  is 


Y 


V8)Ve 


(7) 


6 


and  with  V§/c2  >  2veitopb(z)  vb/cw|e  the  peak  growth  rate  is 


y,  =  2v  .u)2,  (z)V2/v2ui2 
4  ei  pb  b  0  pe 


(8) 


at  k  %  Y4/Vb.  Here  Vq  is  a  measure  of  the  thermal  spread  in  the  azimuthal 
velocity  of  the  beam  ions.  If  krb  <  1,  peak  growth  rate  is  reduced  to 
[l- (k0rfa) _2] y4  because  of  geometry  constraints  and  if  kQrb  <  1,  no  growth  is 
possible.  The  critical  wavenumber  k0  =  /2oipb (z) v^/cVq [see  Eq.  (35) ,  Ref.  18]. 

If  there  is  wave  growth  at  a  given  point  in  z,  the  number  of  e-folds  that 
occurs  is  just 

(V*> 

N  =  \  D  Y(t)dt  ,  (9) 

Jo 

where  the  time  dependence  of  the  plasma  parameters  must  be  considered.  Most 
importantly  vei  decreases  at  Te"3/2(t)  as  the  plasma  is  heated  by  the  passing 
beam.  Note  that  beam  thermal  effects  do  not  completely  quench  the  instability 
because  of  the  finite  plasma  resistivity.  Given  information  from  MHD  consid¬ 
erations  on  the  time  variation  of  the  plasma  parameters  due  to  beam  passage, 
one  can  then  determine  from  Eq.  (9)  whether  significant  azimuthal  beam  current 
bunching  occurs. 


3.2b.  Beam-Weibel  Mode 


For  the  beam-Weibel  instability  the  wave  vector  is  in  the  radial 
direction.  In  this  case  the  spread  in  perpendicular  velocities  associated 
with  the  betatron  motion  of  the  beam  ions  reduces  the  growth  rate  of  the 
instability.  A  measure  of  the  radial  velocity  spread  is  given  by  Vqo^,  so 
that  the  peak  growth  rate  is  given  by 


Y  =  24/5.  .vu  , 
5  ei  pb 


v  (z)/u!2  a8/5 


pe  m 


(10) 


7/5 

at  k  %  y  5/vbam  •  Again  the  number  of  e-folds  that  occur  at  a  given  point 
in  z  can  be  determined  by  Eq.  (9)  and  MHD  considerations.  Since  the  thermal 
spread  in  the  radial  velocity  is  typically  comparable  or  larger  than  the 
thermal  spread  in  the  azimuthal  velocity,  the  beam-Weibel  instability  is 
generally  less  dangerous  than  the  beam-whistler  mode.  A  superposition  of 
both  unstable  modes,  however,  will  lead  to  beam  f ilamentation . 


3.2c.  Plasma-Electron  Mode 


If  beam  thermal  effects  or  geometry  constraints  prevent  the  beam 
from  driving  strong  EM  wave  growth,  instability  can  still  result  from  the 
electron  return  current  established  with  the  plasma.  Because  of  the  high 
collisionality  (vei  >  wce)  of  the  plasma  electrons,  radial  (k  =  ker)  and 
azimuthal  (k  =  keg)  EM  modes  driven  by  the  plasma  electrons  are  essentially 
indistinguishable  (aside  from  small  geometric  effects) .  The  analysis  must 
include  electron  thermal  effects,  as  well  as  collisional  effects  since  the 
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electron  streaming  velocity  is  subthermal  even  early  in  the  pulse  before  beam 
heating  occurs. 

At  early  times  in  the  pulse  when  n>peue/cvei  >  1,  the  peak  growth  rate  is 

y  =  (w  .V  /c)  (2-2-1)1*  .  (11) 

pi  e  p 

Later  in  the  pulse  after  the  beam  heats  the  plasma  if  wpeue/cvei  >  1,  then 
thermal  effects  slows  the  growth  of  the  mode.  The  maximum  growth  rate  is  then 
given  by 


Y  =  2_3/2gj  V2/cu  .  (12) 

pe  e  e 

At  all  times  in  the  pulse  geometry  constraints  prevent  wave  growth  if 
k^b  <  1  where  in  this  case  kQ  =  UpeVe//2uec .  Since  again  this  mode  does  not 
convect,  the  number  of  e-folds  at  a  given  point  in  z  can  be  obtained  from 
Eqs.  (11) ,  (12)  and  (9)  and  MHD  considerations. 

3. 2d.  EM  Stability  Conditions 

If  the  growth  rates  of  the  EM  instabilities  are  slow  enough  such 
that  less  than  one  e-fold  (N<1)  occurs  during  beam  passage,  wave  growth  will 
be  washed  out  by  the  MHD  changes  in  the  beam-plasma  system.  Only  wave  growth 
with  N>1  will  affect  beam  transport.  In  order  to  prevent  the  beam  from 
filamenting,  the  beam  should  then  have 

n.  (z)  rTb(z)/  Z  X  .  \ 

a2  >  V2/V2  >  1.6X1CT9  — -  I  (  )  dt  ,  (13) 

m  *  b  yb  J  \T  V2/ 

o  e 

where  again  Te  is  in  eV's  and  n^  is  measured  in  cm-3.  This  condition  is 
derived  from  Eq.  (8)  for  azimuthal  current  bunching  which  is  typically  more 
severe  than  the  condition  for  radial  current  bunching  since  (%,  is  usually 
greater  than  Vg/V^.  Here  it  is  also  assumed  that  t/2u)pfc (Z) r^Vb/cVg  >  1,  which 
is  the  usual  case.  In  order  to  evaluate  Eq.  (13) ,  a  knowledge  of  the  MHD 
response  to  beam  passage  is  required,  in  particular,  one  must  know  the  time 
history  of  the  electron  temperature. 

If  the  beam  is  warm  enough  to  prevent  beam  f ilamentation  (i.e., 

Eq.  (13)  is  satisfied),  it  is  still  possible  for  the  channel  to  filament. 
Channel  f ilamentation  is  avoided  if 


(14) 


Here  Up  =  mp/mH,  Jb(z)=Z]3en]:)(z)Vj;>,Jj:,  is  measured  in  sA/cm2  and  t. 
is  defined  as  the  time  it  takes  for  the  beam  to  heat  the  plasma 
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to  Tj  =  l^xKT^Z^h in  eV.  If  Te(0)  >t|,  then  tj  =  0  and  the  first 
term  in  Eq.  (14)  does  not  contribute.  If  t^  S  Tb,  then  the  second  term  in 
Eq.  (14)  does  not  contribute  and  the  limits  on  the  integration  of  the  first 
term  go  from  0  to  Tb(Z).  In  Eq.  (14)  it  was  assumed  the  uuperbVe/2uec  >1 , 
which  is  the  usual  case.  Eq.  (14)  sets  an  upper  limit  on  transportable  beam 
current  density  above  which  channel  f ilamentation  will  develop. 

4.  CONCLUSIONS 

The  stability  conditions  in  Eqs.  (5),  (13)  and  (14),  identify  important 
constraints  on  beam  propagation.  Eq.  (5)  together  with  Eq.  (6)  states  that 
the  plasma  cannot  be  heated  above  a  certain  critical  temperature  without 
generating  electrostatic  microturbulence.  Such  turbulence  will  degrade  beam 
quality  and  confinement.  In  general  the  most  severe  constraint  on  Te  exists 
at  the  tail  of  the  beam  before  bunching  occurs  where  the  critical  temperature 
is  defined  by 


Tc(eV) 

e 


H  10 


-6 


p2n3Z5X2 .a8 
b  p  p  ei  m 

Zbnb(z=0) 


1/3 


(15) 


where  np  is  given  in  Eq.  (6) .  By  increasing  np  the  critical  temperature  can 
be  raised  at  the  expense  of  increasing  the  beam  energy  loss  during  transport. 
Operating  at  twice  the  density  specified  in  Eq.  (6)  increases  T^  by  a  factor 
of  two  but  only  increases  the  beam  energy  loss  by  25%.  The  actual  temperature 
that  the  electrons  will  reach  for  a  given  beam  and  background  gas  is  not 
always  easily  estimated  and  must  be  calculated  on  a  case  to  case  basis  and 
compared  with  the  critical  temperature  defined  in  Eq.  (15) .  Here  am  is  a 
measure  of  the  spread  in  axial  velocities  in  the  beam,  which  for  a  monoener- 
getic  beam  is  given  by  AVz/Vb  'v-  2.  If  the  beam  has  a  spread  in  energy  such 

that  <  AE/E,  then  in  Eq.  (15)  should  be  replaced  by  (AE/E)4. 

Growth  of  electromagnetic  waves  can  be  held  to  a  tolerable  level  by 
allowing  for  a  reasonable  spread  in  beam  perpendicular  energy  and  limiting 
the  beam  current  density.  If  the  beam  is  too  cold,  beam  f ilamentation  can 
occur  during  beam  transport.  The  condition  for  good  beam  transport  is  given 
in  Eq.  (13)  .  If  the  beam  does  not  filament  but  the  beam  current  density  is 
too  high,  the  return  current  in  the  channel  can  cause  channel  f ilamentation . 

In  order  to  prevent  this  from  occurring  Eq.  (14)  should  be  satisfied.  The 
conditions  given  in  Eqs.  (13)  and  (14)  may  actually  be  too  severe  since 
channel  hydrodynamic  effects  have  not  been  included  self  consistently  in  the 
analysis.  Channel  expansion  due  to  JXB  force  may  tend  to  prevent  growth  of 
these  transverse  modes.  These  effects  are  presently  under  investigation. 
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IDEAL  DIODE  VOLTAGE  WAVEFORM: 
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<H t)=  - Zt—z  i  0<  t  S  rb<t 

d-Vta)2 


Fig.  1 -Schematic  of  ion  diode,  focusing  region  and  transport  channel  for  light  ion  beam  ICF  module. 
The  ideal  diode  voltage  waveform  for  axial  bunching  is  also  shown  where  <f>0  in  the  initial  diode  voltage, 
ta  is  the  beam  arrival  time  at  the  target  and  fb  is  the  beam  pulse  duration  at  injection. 
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Fig.  2— Schematic  of  beam  propagation  in  a  /-discharge  channel.  Here  'tnijVj,  -  Zde0  (where  Zd  is  beam 
ion  charge  state  in  the  diode  and  o  is  the  diode  voltage),  am  i>,  the  maximum  injection  angle,  r5  is  the 
beam  spot  size  (or  aperture  si/.e)  and  rh  is  the  beam  radius  in  the  channel.  Also  shown  is  the  channel 
current  required  to  confine  the  beam  where  7.h  is  the  charge  state  ot  the  beam  ions  in  the  channel  and 
Fb  = 


PLASMA: 

Te&Tj 

^ei  ( Te )  >  ^ce 

J/2 

Ve“ZbnbVb/ZPnP<  (T%ne) 


BEAM: 
Vb/c  « I 
nb«  np 
Vt 

—  w  a 

Vb 


Fig.  3-Schematic  of  particle  distribution  in  vz.  Here  i,  e  and  b  are  the  plasma  ions,  plasma  electrons 
and  beam  ions,  respectively,  and  V,  is  the  measure  of  the  spread  in  the  beam  distribution.  Other 
typical  beam-plasma  conditions  are  also  listed. 
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